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' Extending our recent work (hep-th/0310106) we study the nonsinglet sector of c = 1 matrix 
model by renormalization group analysis for a gauged matrix quantum mechanics on circle with 
^ [ an appropriate gauge breaking term to incorporate the effect of world-sheet vortices. The flow 
equations indicate BKT phase transition around the self-dual radius and the nontrivial fixed 
points of the flow exhibit black hole like phases for a range of temperatures beyond the self-dual 
point. One class of fixed point interpolate between c = 1 for i? > 1 and c = as i? via 
black hole phase that emerges after the phase transition. The other two classes of nontrivial 
fixed points also develop black hole like behavior beyond R = 1. From a thermodynamic study 
of the free energy obtained from the Callan-Symanzik equations we show that all these unstable 
phases do have negative specific heat. The thermodynamic quantities indicate that the system 
does undergo a first order phase transition near the Hagedorn temperature, around which the 
new phase is formed, and exhibits one loop finite energy correction to the Hagedorn density 
of states. The flow equations also suggest a deformation of the target space geometry through 
a running of the compactification radius where the scale is given by the dilaton. Remarkably 
there is a regime where cyclic flow is observed. 
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1 Introduction 



The c = 1 matrix model have been proved to be very powerful in describing the two-dimensional 
string theory to all genus, both in string perturbation theory and in a nonperturbative sense 
Q 121 ini HI- The underlying rich structure of the non-relativistic quantum mechanics of free 
fermions (exhibited by the singlet sector) makes the theory solvable to a very high degree and 
virtually renders any quantity calculable to all orders in genus expansion. Recently the c = 1 
matrix model is realized as the effective dynamics of D-branes in c = 1 non-critical string theory 
which exhibits the duality between c = 1 matrix model and 2D quantum gravity coupled to 
c = 1 matter as an exact open/closed string duality jSlinilZllHliniliniinillSllIHl. The quantum 
mechanics of SU{N) invariant matrix variables in an inverted oscillator potential is visualized as 
the quantum mechanics of open string tachyons attached to unstable DO branes that decays 
into {i.e. dual to) Liouville theory coupled to c = 1 matter describing 2D closed string theory 
together with its DO branes. 

The random surfaces mapped by the c = 1 matrix model can also be embedded in a circle 
of radius i? as a compactified Euclidean theory or equivalently as a Minkowski signature string 
theory at finite temperature, where the free fermion representation is not sufficient due to the 
active role of the angular degrees of freedom belonging to the nontrivial representation of SU{N) 
|14 | I15 | [T6]. From matrix quantum mechanics analysis the states in the nontrivial representation 
of SU{N), the nonsinglet sector, are understood to correspond to vortices on the world-sheet 
with wave functions given by the Young Tableaux of SU{N). The number of boxes counts 
the vortex charge. Restricting the theory only to the SU{N) singlet sector gives rise to the 
continuum limit where the effect of the world-sheet vortices are absent and the sum over the 
random surfaces obey T-duality Pl]. Keeping the vortices triggers Kosterlitz-Thouless phase 
transition at self-dual radius that breaks the R 1/R duality symmetry. Thus the states in 
the nontrivial representation of SU{N) are important in understanding the dynamics of the 
world-sheet vortices and are expected to give rise to interesting phenomena like formation of 
2D black hole [HI. 

An interesting solution of the two-dimensional string theory, apart from the flat space 
background with linear dilaton, is the two-dimensional black hole JSl ^1 1201 ■ Some of the 
early attempts to get matrix-model description of two-dimensional black hole are described in 
ini 1221 123 EH 123 ESI EZl EH- it was partially believed that the nonperturbative formulation 
of two-dimensional string theory in terms of an integrable theory of noninteracting nonrela- 
tivistic fermion representation of the matrix quantum mechanics can deal issues like black hole 
evaporation and gravitational collapse. But it turns out to be remarkably difficult to find a 
matrix model description of the two-dimensional black hole due to lack of clear understanding 
of the target space physics of linear dilaton background in the c = 1 matrix model. The wave 
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equation should carry information about the black hole supplemented by the dilaton. In free 
fermion representation, the time- of- flight coordinate r, related to the matrix model coordinate A 
by A = y/ 2 fia' cosh T, cannot be identified with the Liouville coordinate although this identifi- 
cation arises naturally from the collective field theory approach [211 EDI- This is due to the fact 
that the eigenvalue coordinate A has no obvious geometrical interpretation on the discretized 
world-sheet. According to the correct identification is through the loop operator W{1), 

which has a clear geometrical meaning on the world-sheet, by an integral transform. Also it has 
been pointed out in that exact solvable structure, especially the W{oo) symmetry of the 
c = 1 model makes black hole hard to describe. 

One would hope to understand the situation better by working in a more general represen- 
tation including the world-sheet vortices. However, the study of world-sheet vortices is hard 
due to lack of solvable structure. In the spherical approximation, they could be studied in 
dual matrix description considering discrete time jHHESlEni- In [Hli using twisted boundary 
condition (j){27rR) = Q''(j){0)Q, the partition function was studied in a given representation for 
the standard matrix oscillator (with a stable quadratic potential). The partition function in 
the presence of the adjoint representation (a vortex-antivortex pair) in the double scaling limit 
was then studied by analytically continuing to the upside down oscillator . However, a direct 
analytical continuation is not possible as the standard oscillator has larger symmetry than the 
upside down one and does not have any information about the cut-off provided by the interac- 
tion terms in the matrix potential. Hence the analytical continuation had to be completed (in 
the spirit of [37]) by a suitable guess about the cut-off dependence and was argued to work at 
least for the adjoint representation. 

Based on the above approach of connecting the world-sheet vortices with the nonsinglet states 
of the matrix quantum mechanics, an integrable system has been constructed which is an inte- 
grable Toda chain hierarchy interpolating between the usual c = 1 string and the Sine-Liouville 
background [IZj^. Using the duality conjecture by Fateev, Zamolodchikov and Zamolodchikov 
(FZZ) [21], it was proposed to be a matrix model (the KKK model) for two-dimensional black 
hole at i? = 3/2 that relates the two-dimensional black hole background {SL{2, R)/U{1) coset 
CFT) to the condensation of vortices (interpreted as the winding modes around the Euclidean 
time). The effect of the winding modes on the world-sheet were incorporated in the matrix 
quantum mechanics path integral by integrating over the twist variables (the matrix holonomy 
factor around the compactification circle) along with an appropriate measure and the twisted 
partition function. The FZZ correspondence allows one to consider the winding mode conden- 
sation from the Sine-Liouville side to construct the appropriate matrix model for the black hole 
background, avoiding dealing directly more complicated black hole geometry. However, not 

^The description is related by T-duality to the Toda integrable structure of the c = 1 string theory perturbed 
by purely tachyon source [SHI- 
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considering the black hole background directly has some disadvantages for the following reason. 
Although FZZ conjecture has been tested by calculating various correlators, it is not clear how 
to get information about the black hole metric from the Sine-Liouville side. On the other hand, 
considering the thermodynamics of two-dimensional string theory above the temperature corre- 
sponding to the Berezinski-Kosterlitz-Thouless (BKT) phase transition |^ , it has been 
proposed J7] that the nonsinglet states (or the vortices on the world-sheet) fills the Hagedorn 
density of black holes p{E) ~ E~'^^e^"^ , since the perturbative spectrum of two-dimensional 
string theory has very few states, namely the massless tachyon. The Hagedorn spectrum of 
states should then be obtained by a direct counting of the nonsinglet states in the c = 1 ma- 
trix model. From a Hamiltonian point of view, to count the states in a given energy interval 
one needs to diagonalize the Hamiltonians for representations corresponding to the large Young 
tableaux. This is difficult due to Calogero type of interaction of the eigenvalues with the SU {N) 
spin structure. Also, there are other important questions like understanding black hole phases 
at any radius other than the fixed radius R = 3/2 (to which the unstable black hole can decay 
to), the possibility of observing Hawking radiation and the like. 

The main motivation of the present paper is to address above questions from an explicit 
study of the nonsinglet sector and the BKT phase transition directly from the matrix quantum 
mechanics path integral with periodic boundary condition. Instead of working in any particular 
representation let us have recourse to the renormalization group approach that we developed in a 
recent paper 13j for the c = 1 matrix model on a circle of radius R. There we described a detail 
analysis of how the two coupling constants in the double scaling limit with critical exponent 
fiow with the change in length scale. The motivation came from an earlier work of Brezin and 
Zinn- Justin [43^ on large RG analysis of c = matrix model. The scheme reproduces the 
known results of the solvable sub-sector of the matrix quantum mechanics, namely a non-trivial 
fixed point with the correct string susceptibility exponent of the c = 1 model, T-dulaity and 
the expected logarithmic scaling violation of the free energy. Also it exhibits, qualitatively, 
the physically interesting situations due to the effects of nonsinglet sector, which can not be 
simplified because of the lack of the solvable structure. For example, from the running of the 
prefactor of the partition function, written in the renormalized couplings, analogous to the 
running due to the wave function renormalization, the free energy is observed to change sign 
near R = 1 for small value of the critical coupling. This is reminiscent of the BKT transition 
at self-dual radius triggered by the liberation of the world-sheet vortices. We would like to 
understand the detail nature of the nontrivial fixed points of the fiow that describes the physics 
beyond this transition. 

To capture the effect of vortices on the fiows and the fixed points more clearly and to 
introduce a new coupling that would act like vortex fugacity, in this paper we analyze the 
behavior of the following gauged matrix model with simple periodic boundary condition and 
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with an appropriate gauge breaking term 



7</,jv(27ri?)=<^jv(0) 

- (iV) Tr ^"^""rft {\{DMt)f + I'Pl^) - + ^} 



where the covariant derivative D is defined with respect to the pure gauge A{t) = Q{tyQ{t) by 
D(j) = dt4> + Unhke taking the course of integration over the twist fields with a proper 

measure incorporating winding modes around the Euchdean time in the path integral (as con- 
sidered in fZj), here the integration over the gauge field A(t) = Q^Q(t) with an appropriate 
measure provided by the gauge breaking term inserts world-sheet vortices in the partition func- 
tion where a acts as the vortex fugacity. This can happen through the insertion of operator of 
the form exp(— aJ^) that counts vortex number, where ~ nN. Without the gauge break- 
ing term the system is projected to the singlet sector. We observe that one class of nontrivial 
fixed points of the fiow give rise to a pair of c = 1 fixed points at large R with one unstable 
direction. As R is decreased the fiow passes through Kosterlitz-Thouless phase transition at 
R = 1.03, where the operators coupled to the vortex fugacity become relevant indicating the 
liberation of the world-sheet vortices as expected at the self dual radius R = 1 ,15^- Between 
0.67 < R < 0.7 the pair of fixed points become purely repulsive fixed points of large coupling 
and exhibit negative specific heat and one loop correction to the Hagedorn density of states very 
similar to those exhibited by an unstable Euclidean black hole in fiat space time. The change 
of entropy exhibits a discontinuity at R = 0.73, little above the BKT temperature, indicating 
the Hagedorn transition to be first order. As R is further decreased the fiow ends up into a pair 
of c = fixed points. The other two classes of fixed points also show black hole like behavior 
beyond the BKT phase transition. This indicates he existence of other black hole phases at 
other radii of compactification. Also a running of the compactification radius with the scale, 
thought as dilaton, suggests a deformation of the target space geometry that might be crucial 
to visualize those black holes. We observe cyclic fiow below the BKT temperature presumably 
due to resonance of high spin states indicating stringy behavior In fact, beyond the phase 
transition, where different nontrivial fixed points exhibit black hole like behavior, the cycles 
become very complicated. The phase structure probed by the RG analysis thus remarkably 
captures the expected properties and looks promising in understanding the much unexplored 
physics of the nonsinglet sector. The dynamics in the neighborhood of the black hole like fixed 
point needs to be studied in detail to understand its explicit nature. In this paper we will 
motivate such studies regarding the role of the nonsinglet sector from the observations from our 
Renormalization group analysis. 
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The plan of the paper is as follows. In section 2 we describe the RG calculation for our 
model (jl-H) . In section 3 we analyze the flow and discuss the phase structure. We observe the 
BKT phase transition near the self dual radius (or a Hagedorn transition little below the BKT 
radius) and existence of a black hole like fixed point past the transition. In section 4 we analyze 
the thermodynamics of the black hole like fixed point and comment on the sense in which it 
resembles the thermodynamics of an unstable (thermal) Euclidean black hole. In section 5 we 
conclude with discussion and open questions. 



2 The Large N RG Calculation 



In this section we will summarize the renormalization group calculation for the c = 1 gauged 
matrix model on a circle, with a gauge breaking term appropriate for capturing the nonsinglet 
physics. The details of the method for the ungauged model with periodic boundary condition 
is studied in [13 . 

Let us consider the partition function for the + 1 dimensional matrix variables 



ZN+i[g,a,R] 
exp [ - (iV + 1) Tr 



</.jv+i(27rR)=</.jv+i(0) 



dt 







a 



The covariant derivative D and the gauge field A{t) are defined respectively as 

D(f) = dt(j) + [A, (f)] , A{t) = n{tyn{t) , n e u{n) . 

Expanding the covariant derivative, the partition function is rewritten as 



ZN+i[g,a,R] 



- (iV + l)Tr p'^ dt {l<PN+i{tr + ^0^+,(t) - |0^+i(t) 

}] 



'0jv+i(27rR)=</,]v+i(O) 

exp 



1 A"^ 

+AN+l{t) [0^+1 (t), 0^+1 (t)] + - [AN+l{t),(PN+l{t)]' + 

2 a 



(2.1) 



(2.2) 



(2.3) 



The A]\r+i{t)[(f)N+i{t), (i>N+i{t)] term above is crucial to study the nonsinglets. Even though 
they are present, the gauge invar iance tries to project the system to the singlet sector while the 
gauge breaking term prevents to do so. In jT3], a finitely large radius representation of singlet 
sector was obtained by throwing this particular term by hand as the nonsinglets are confined 



7 



at small temperature. In JB], the partition function for one vortex/anti- vortex pair, i.e. in 
the adjoint representation was calculated by analytical continuation from the twisted partition 
function of the standard harmonic oscillator to that of the upside down oscillator. For a = 0, 
the gauge fields are forced to vanish and the partition function reduces to that of ungauged 
matrix quantum mechanics on circle. 

Because of the gauge breaking term, the integration over all possible configurations of A{t) 
formally inserts (the gauge invariant) operator expTr(— aJ^) in the partition function. 



dA expTr^ \-2iAJj ~ exp(— aJ^) ~ exp(— A^an) 



(2.4) 



Here is proportional to the quadratic Casimir invariant C{n) ~ Nn. Characterizing the 
irreducible representations in terms of the number of the white boxes n in the Young tableaux, 
the quadratic Casimir only depends on n to the leading order in A^. The reason for this behavior 
is that, expTr(— aJ^) acts on the states |Adj) in the adjoint representation (belonging to the 
nonsinglet sector) of the MQM in the gauge invariant way, 

A- J|Adj)" = a n|Adj)". (2.5) 

The parameter a behaves like the fugacity of vortices. The operator expTr(— aJ^) therefore 
counts the vortex number. 

2.1 Integrating out a column and a row of the matrices: 

Now we decompose the (A^+ 1) x (A^ + 1) matrices into N x N blocks and A^- vectors and scalars 
as follows 

, /,x ( (t>N{t) VN{t) , 

(t)N+iit) , (2.6) 

and 

A-^.M = i^^fl . (2.7) 

The scalars e and t] are of relative order 1/A^ and can be ignored in the double scaling limit. Even 
though the matrix elements are not independent (0 being hermitian and A being antihermitian) 
one can always choose such decomposition for each value of A^ which does not prevent the action 
to be written in the same form in terms of the covariant derivatives. The resulting partition 
function can be written as 



8 



exp " " " ' ^ 

'0 



- (iV + 1) Tr rft + -<P%{t) - ^-<P%{t) 

+A^{t) [Mt),Mt)] + I [AN{t),Mt)? + ^}] ■ 
X [ P^a(^)P^a*(^)P^^;(^)P^^;*(^) 

Jv,v*(2-itR)=v,v*{0) 

/2-kR 
dt \^v*v + v*Anv - v*Anv + v*{A% + 1 



9(t>N)v 



a 



+a* ( H — )q + i^*(0Af + 2/ljv0Ar — 0Ar^Ar)a — v*(l)Na 



+a*{—(j)N + '2(I)nAn — AN(j)N)v + a*(t)NV 
-\-a*va*v + v*av*a — v*va*a — a*vv*a^ 



(2.8) 



vit) 



Rescaling the vectors v(t) —>■ , 

^ ^ ' ■sJ2-kR(N-^\) 

the partition function turns out to be 



and a(i) 



a(t) 



^2i^R{N+\) 



, the V and a dependent part of 



I[g, a, 4>N, An, R] 



V''a{t)V''a*{t)V''v{t)V'\*{t) 



'v,v*{2nR)=v,v*{0) 

- J ^ [v*v + v*Anv - v*Anv + v*{Al + 1 - gcf)N)v 



exp 



+a*{-4>N + 20ArAAr - AN(f)N)v + a*(f)NV 



27rR{N + 1) 
Using the foUowing Fourier transformation 



— r(a va V + V av a — V va a — a vv a) > \ . 



(2.9) 



m=—oo ^ 



with 



dt , (n-ra) 

e 



,0 27ri? 

the (v, a) -integration can be expressed as 



, 5(t - 1') = ^ y , 



(2.10) 
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I[g,a,(j)N,AN,R] 



{^2tiR{N + 1)) 



AN 



n ^"'*mda*m) ( n 



exp 



exp 



exp 



exp 



exp 



E 



m,/ k 

Yl ^*mCZl (t>m-l-k<Pk)ai } 
m,l k 

/ (m — I) m \ ^ , 



m—l—kYk 



E f {m — I) m X ^ X ^ 



m,l 



(2.11) 



where we have neglected the 0{1/N) terms. The above integration, as an extra piece to the 
matrix integral, contains the vectors v{t) (or quarks) generating boundaries on the Feynman 
diagrams. In fact in vanishing a limit, this is essentially equivalent to the model considered in 
|45| ITH] . For generic value of a, the above integral should insert nonsinglet boundaries on 
the world-sheet that arise after the BKT phase transition. 



2.2 One loop Feynman diagrams 

In order to carry out the v and a integration diagrammatically, let us now define the following 
operators 



mn 

Ot-M.A) 



(5r, 



a 



[m + l) 



R 



Am-l + ^ ] Afn-l-kAk 



>m-l-k<i>k 



.{2m -I) 

■ 1 

R 

{21 - m) 



^m—l 



^m—l 



+ 2 X Am-i-k(l)k - X 

k k 

+ 2 X^ (pkAm^i^k — X^ Am-l~k(Pk 

k 



m—l—k 



[2.12) 



The inverse of these operators define various propagators and vertices according to figure ^ 
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v*v -1 ^ a*a -1 

[O ] = " [O ] 

L mm J L mm J 



a"^ a 



mm-' m jn mm- ^ ^ 



v*v -1 

[0„,] - 



m 



a*a -1 



m 



v*a -1 

[0„_,] - 



m 



a*v -1 

[0„,_,] - 



a"^ 
m 



Figure 1: The propagators and vertices. 



Hence the integral becomes 



I[g,a,^r,,Ar,,R,N] = ^-^—^ exp [ - { ^ (0:;7„ 0, A))-^ 

mi .li 



"12,^2 "13,^3 



where the gaussian part is as follows 



Performing the gaussian integration, we get 



(2.13) 



Jo[a, i?, iV] = 1 ( n ( n ^^1^^^) [ - + <^'^nym)] ■ (2.14) 



/o[a, R, N] = (27ri?)2^vr2^ exp [ " Trj ^ In + l) 1 

m 

Inserting this into (j2.13|) . the (f , a)-integration becomes 



exp 



Tr<i ln(-)l 



I[g, a, (pN, An, R, N] = C{a, R, N) T,[g, a, (pN, An, R, N] 



where 



C{a,R, N) 



ir^Ra 



(A^+ l)2sinh7ri? 
11 



N 



UJrn/Rr^ 



(2.15) 



(2.16) 



:2.17) 



In ()2.16|) . T,[g,a, (f>N, An, R, N] represents sum of one loop Feynman diagrams as shown in figure 

121 






2.3 Evaluation of the diagrams 

We evaluate the diagrams to express T,[g, a, (pN, A^, R, N] as a series expansion in the couplings, 
g, a, and the Fourier modes of the fields 0^, ^m- The full series expansion is given in the 
Appendix A. We will now discuss how to express this series into an expansion in g, a, $(t) and 
A{t) by suitable inverse Fourier transform. 

In order to evaluate the one loop correction to the effective action, we inverse transform 
the Fourier modes according to the rule ()2.im) and sum-up the set of infinite series using the 
formulae discussed below. 



E°° exp[i{m/R)t\ ttR cosh.{TT aR — at) 
5 = ; 7^ , < t < ZttR , 
^ + a sinh iraR ' - - 

m=—oo 

m/R exp[z(m/i?)t] mR sinh(7rai? — at) 

^ + smh T[aR 

m=— cxD 

inR , t — > , 

—mR , t 2nR , 

0, t = = 27iR. (2.18) 

Because of the discontinuity in the last series shown above, all other series, which look like 
higher derivatives of the above functions, has delta function or derivative of delta function like 
behavior. This behavior actually dominate over the regular part of those functions and hence 
dominate the contribution when the functions are integrated over : 
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f rnVR^eM^{^/R)t] ^ ^^^^^^^ ^ 0<t<2nR, 



m=— oo 



oo 



v-^ rrr/R^ cxp[i(m/i?)tl n „ ^ 

2_, — 2 -27riRS'{t) , 0<t< 2'kR, 

m=—oo ^ ^ 

£ "^'/"^'jfP^il;"/-^^^] ^ -2nR5"(t), 0<t<2nR. (2.19) 

m=— oo i?2 ~l~ ^ 

One way to see the above behavior is as follows: 



m=—oo 
oo 



E Pfr. - E 1 = 2.i?<5(0) . (2.20) 



-52 + Oj 
m=— oo it^ 



Also 



"55 T O 

m=— OO ii m=— oo 

oo Q / T->o oo 

m=— oo /[^ m=— oo 

|, mVfi!j«M . E |jexpM*l = -2.firW, 0<t<2.fl. 

m=— OO m=— oo 

V - E = -2nR5"{^) . (2.21) 

m=— oo it'' m=— oo 

Note that while integrating over these functions, we have to take into account their periodic 
nature and the intervals over which they are defined. We also have to deal with subtleties when 
these intervals are strictly inequalities. 

These sums give the corrections to the coefficients of the various terms in the action, after 
T\g, a, 0JV, Am-i R, N] is exponentiated and log-expanded around = (small field approxima- 
tion). Since after doing the inverse Fourier transform, the various terms containing (f){t) and 
A{t) has nonlocal integrals over several one dimensional dummy time variables, we breakup the 
variables into center of mass and relative coordinates. Then we expand the functions about the 
center of mass coordinates, assuming the relative coordinates to be small enough, and consider 
integration over the relative coordinates. 
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Now taking into account all these considerations, after evaluating the integrations over rel- 
ative time variables (see Appendix B), the expression for T^lg, a, (pN, ^at, R, N] becomes 



E[g, a, cPn, An, R,N] = 1 + Fgi{R) g / dt Tr0(t) 

+{Fg2{R) / + F^2iR) a} I dt Tr-02(t) + {p^^iR) g^ + F^^m «} / dt Ti-4>\t) 

+G2iR) / dt TrA^{t) + {GgsiR) g^ + Ga^{R)a} / dt Tr^(t)0(t)0(t) 
Jo Jo 

+{G',s{R) g' + GURh} / dt Tr A{t)^{t)<f>{t) + {F,s{R) g' 

Jo 

p2i^R -1 p2-KR 

+Fags{R) ag} / dt Tr-^^t) + {G^^iR) « + Gg^{R) / dt Tr A{t)(l>{t)A{t)(l>{t) 
Jo "J ^0 

+{GUR) « + G'g,{R) g'} / dt Tr A'{t)<f>'{t) . 

Jo 

(2.22) 

The functions F{R)s and G{R)s are defined as follows 



Fgi{R) = cothTri?, Fg2(i?) = ^(1 + cothVi?) + ^coth7ri?cosh27ri? 

Fa2{R) = -l/(^/?) + 8-4coth7ri?, F„2(i?) = -4 cothTri? , 

Fg2(R) = \ (1 + 87r^i?2) sinh AnR - cosh 2t^R - — cosh AnR) , 

sinh vri? 64 6 16 

ttR ttR 

G2(R) = -3coth7ri?+ ^ — (1 + cothVi?) , 

2 sinh ttR 2 

ttR 

F„3(R) = ^ [47ri?(3cosh7ri? + 2cosh37ri? + 2cosh57ri? 

^ 64sinh^7ri?(cosh27ri? + cosh47ri?)^ 

+ cosh 7t:R) + sinh ttR + sinh BttR + 2 sinh TttR + sinh dnR] , 

„ , 397ri? , , 9 „v 33 , „ 

^^3(7?) = G'g3(i?) = —{ttR cosh 27ri? - - sinh 27ri?) , 

Ggi{R) — ('k'^R^ cosh 37ri? + sinh 57ri? + sinh 37ri? — \ sinh 7ri?) , 

8 sinh 7ri? 8 4 8 

Gaz{R) = G';,3(i?) = -^C0th7ri?+— -^(^ cosh 7ri?+-^ sinh TTi?), 

3 smh Trit b 12 

^^^(i?) = -2coth7ri?, G';,4(/?) = -cothTTi?, G^4(i?)=0. (2.23) 
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2.4 Elimination of the tadpole term 

The term proportional to J^""^ dt is unwanted. To remove this term, we change the back- 
ground by (j){t) (f){t) + /, and set the net coefficient of the hnear term to zero. This fixes the 
value of / to be 



- + N ) iV + N N ) 

As /+ ~ 0(1) and /_ ~ 0(1/ N), we choose the smaller shift 

/_ ~ gFg^/N (2.25) 

to suppress the contribution from the higher order terms. Note that the coefficient of 0"'(t) 
term would contribute a term proportional to / in the coefficient of 0"~^(t) term. Thus the 
coefficient of (6'^(t) in the coupling would have an 0{1/N) contribution in the coupling 
from the (f)'^{t) term, if we had turned it on. Also the contribution from the [A(t), 0(t)]^ term in 
the coefficient of the A'^{t) term can be ignored as it is of 0{1/N'^). After accommodating all 
such changes, the expression for Z^+i becomes 



//■ 
V(f)VA exp -NTi dt{{l + l/N 

-{aF^2 + g'Fg2)/N)^(t)/2 + (1 + 1/N - (aF^^ + Q^Fg^ + Fg^)) / N)(t>\t) /2 
-{g + g/N + {agF^gs + g'Fgs)/N)<p%t)/3 + (1/a + l/{aN) - G2/N)A\t) 

+(1 + 1/iV - (aG^g + g^Gg3)/N)A{f)mm 
-(1 + 1/N + {aG'^3 + g^G'g3)/N)A{mmt) 
+(1 + l/N - (aGa4 + g''Gg4)lN)A{t)(t>{t)A{t)(t>{t) 

-(1 + 1/A^ + {aGU + g'G'g,)/N)A\t)ct>\t)]\ . (2.26) 
Here the expression for =^(q;, g, R, N) is given by 

- g^F^jN^ + 0{1/N^) . (2.27) 
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2.5 Rescaling of the fields and the variables 

We will now rescale the fields and the variables (time t and the conjugate momentum 1/R) to 
restore the original cut-off (in the Wilsonian sense): 

t t'{l + hdl), R-^ R'{1 + hdl), 
<P{t) ^ pcP\t') , 

A{t) ^{l-h dl) r]A'{t') , (2.28) 

where, 

dl = 1/N , h = h{R) + ^ Cij g'a^hij{R) . (2.29) 

The exact functional form of h can be guessed from the behavior of the Feynman diagrams. We 
will see that h{R) is the scaling dimension of the operator coupled with mass parameter, the 
coefficient of the 0(t)^/2 term (1/a', we have set a' = 1 here), and is appearing in the universal 
term of the beta function equation of the mass parameter. It is interesting that we will also see 
h{R) to appear in combination with other numbers as the scaling dimensions of the operators 
coupled with the cosmological constant g, and the fugacity a. Therefore, being in the universal 
term of the beta function equations of the couplings g and a, h{R) determines the radius at 
which the corresponding operators become relevant and could trigger a phase transition. The 
latter will happen if there is any discontinuous change in the free energy like the flipping of 
sign. In the world-sheet picture one expects a phase transition at the self dual radius R = 1 
due to the liberation of the world-sheet vortices The world-sheet free energy changes sign 
at that radius due to a contest between the entropy of the liberated vortices and the energy 
of the system. We will come back to this in the discussion of the beta function equation of 
the fugacity parameter a and will see that its universal term do reflect such a transition. Also 
through the rescaling relation of R ()2.28j) . h{R) determines the change of the radius with the 
scale and hence will play a role in discussing the thermal properties of the fixed points of the 
RG flow. 

Now in doing the rescaling of time, as shown in ()2.28p . the field A{t) automatically peaks up 
a factor of {1 — h dl) because of the definition ()2.2j) . Since there is no other constraint on A{t), 
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we are here free to rescale it by an arbitray factor 77. The rescaled action therefore looks hke, 

S' = NTt / dt' (1 - hdl)p^(-^'{t')^ + C3 vA'{t')<^'it')^'it') - c'^ 7]A'{t')^'{t')<^'{t') 
Jo ^2 

+C4 r]^A'{t')(!>\t')A'{t')(!>'{t') - r]^A:\t')^'\t'] 

2 

(1 + hdl)p^^<!>'\t') - 1(1 + hdl)p^C3<l>'^{t') + (1 - hdl)c^ "^A'it'f 
2 3 a . 

(2.30) 

The coefficients q and q, z = 2, 3, 4 can be read off by comparing with that of the renormahzed 
action in ()2.2(jj) . Setting the coefficient of the kinetic term 0^/2 to one, p can be fixed as 



p = l + ^[h-l + g^Fg2 + aF^2]dl + 0{df) 



(2.31) 



Similarly setting the coefficients of A{t)[(l){t) , (pit)] and that of the [A(t), 0(t)]^/2 term respec- 
tively to one, we have 

—f] = -rf] = 1 , —f] = —V = 1 • (2.32) 

C2 C2 C2 C2 

In other words, this fixes f] to 



r/ = ^ = 1 + [(G,3 - ^,2)^?' + {Ga3 - F^,)a]dl + 0(t^/2) 
C3 



(2.33) 



along with the constraints 

C3 = C3 , C4 = C4 , C3 = C2 C4 . (2.34) 
Accordingly the coefficients of 0^(t)/2, g(f)^(t)/3, and (l/a)A^(t) respectively are modified as 

C2 + (5c2 = (1 + 2/i c//)? = 1 + [2/i + {Fg2 - Fg2 - Fgi)g^ + (F„2 - Fg2)a]dl + 0{df) , (2.35) 

C2 



C3 



C3 + 5c3 = (1 + 5/1/2 rf/)^ = 1 + [5/1/2 - 1/2 + (3F,2/2 + Fgs)g' + (3F„2/2 + F^gs)a]dl 



+0{df) , 



(2.36) 



and, 



C2 C4 



C4+5c4 = (l+hdl)^ = l + [{l-h)+2{Gg^-Fg2)g'+{2Ga3-2F^2-G2)a]dl+0{dl'). (2.37) 
C3 

Setting the coefficient of (/)(t)^/2 to one (this implies keeping the mass parameter at the fixed 
point with unit magnitude) results in an extra constraint 

^2 = . (2.38) 
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2.6 Beta function equations 

The effective action is of the same form as the bare one, but with renormahzed strength of the 
couphng and the fugacity. The resulting partition function is given by 



ZN+i[g' , a' , R'] = A 



exp 







a' / 



(2.39) 



where 



A'^' = C{R, N) exp[-27rRN^^{g, N, R)]p^^ . (2.40) 

Neglecting the 0{dP) terms, the renormahzed fugacities (couplings) and the vacuum energy 
(the prefactor of the partition function) are expressed in terms of the bare quantities as follows: 



g' = ^+[(5/i/2-l/2)^ + (3F,2/2 + F,3)^=^ + (3F,2/2 + F,,3)«^]rfi, 



= - + [(l-/i)- + 2(G,3 
a' a a 

TT^R^a 



Fg2)^ + {2Ga3-2F^2-G2)]dl 



X' 



1 + ln 



sinh 71 R. 



dl + ^ [{h - 1) + g^Fg2 + aFa2\ dl 



(2.41) 



Here, in simplifying the part C(i?, N)^ in the expression for A', we have assumed that for any 
value of i?. 

In I I — 

.N 



C{R,Ny 



— exp 



sinh ttR 

3 7?2 



.{N + l){n\)^ 
N KsinhnRMJ ^ ^ ' ' 



n— ►oo, W— ►oo 



(2.42) 



Also, the term ex.Y>[-2'KRN'^^{g,N,R)\^ contributes only a factor of 1 as ^{g,N,R) 
0{dP). Hence the resulting beta function equations are expressed as 



Pa 



(1(1/ 

— = -(1 - h)a - (2G'«3 - 2F,2 - G2)a^ - 2g^a{Gg^ - F^^) , 



dX , 



r Ti^R'^a 



sinh ttR 
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(2.43) 



along with the constraints, 



Oi{Fa2 + Ga4, — SGas) + 5'^(-^g2 + Gg^ — 2G g^) 

2h + g\Fg2 - Fg2 - Fgi) + «(F„2 - F^2) 



0, 



0, 



0. 



(2.44) 



The relation ()2.28|) indicates in some sense a running of the radius R 




(2.45) 



This suggests a deformation of the target space geometry if one considers the scale to be dilaton. 
In the above equations, all the functions Fs and Gs are hyperbolic nonlinear functions of R, as 
defined in flT^ . 

Now before going to the detail analysis of the fixed points let us try to understand few 
things about the structure of the beta function equations. The much of the structure depends 
on understanding the quantity h. One can clearly see ioi g = and a = 0, the gaussian model is 
never expected to flow and hence for such a fixed point h = {i.e. corresponding to the trivial 
rescaling t' = t and R' = R). The situation is different for a non- vanishing h. Demanding the 
mass parameter (the coefficient of the (fP'{t)/2 term) to be set at some fixed value (we will 
use 1 for simplicity) and not running, one can easily determine some h = h{R) for nontrivial 
fixed points from the set of the beta function equations and the set of constraints. This will be 
discussed in the next section. This is extremely interesting as it could indicate a phase transition 
at certain radius due to turning on the different operators coupled to the relevant couplings. 
Note that in ()2.43|) the linear term in a in changes sign at h = 1. This indicates that the 
coupling corresponding to the vortex fugacity becomes relevant in the h > 1 region. In the next 
section, we will study this as an indication of the liberation of the world-sheet vortices due to 
Kosterlitz-Thouless type of phase transition undergone by the c = 1 matrix quantum mechanics 
with nonsinglet sector. 

Keeping = 1 for simplicity is consistent with the value of the mass parameter (M^ = ^) 
one originally works with in the matrix partition function to visualize the matrix path integral 
as the generator of the discretized version of the Polyakov path integral of 2D bosonic string. 
In recent identification of the matrix quantum mechanics with the quantum mechanics of open 
string tachyon on unstable DO-branes the mass parameter = ^ is identified with the open 
string tachyon mass. A framework of the fiow of a general {i.e. h = ) could be interesting 
to discuss the presence of the boundaries. 
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3 Analysis of Flow Equations and Phase Transition 



3.1 The fixed points 

The fixed points of tlie flow are given by the simultaneous solution of the beta function equations 
Pg = Pa = . The Gaussian fixed point Aq is given by 

g* = 0, a* = 0, h = 0. (3.1) 

The nontrivial fixed points of the flow AKg*"^ ^ 0, 0), A^(0, a* 0) and A^ig*"^ 7^ 0, a* ^ 0) are 
as follows, 

= ^2/ {Fg2 — Fg2 — Fgl — 2 Fg^ — 3Gg4), ] Ri^ , 

h = —. (^^^ - - Fgl) ^ ^3 2) 

{Fg2 - Fg2 - Fgl - 3G g4 " 2 Fg3) 



where Ri is given by, 
Similarly, 



Gg3 — G'g^ , Gg4 — G'g^ . (3.3) 



^2 = 2/ {Fa2 — Fa2 — SGai + 26*2 — 2Fag3)] -R2 j , 

h = —. , (3.4) 

{Fa2 — Fa2 — SG^i + 26*2 ^ 2^0,^3) 



where R2 is given by. 



Ga3 — , Ga4 — G'^^ . (3-5) 

And lastly. A* (51*2 ^ 0, a* 7^ 0) is given by, 

g*'^ = 2{Fa2 + G'^^ — 36*0,3 + 

I {(-^a2 + G'^^ — 36*03 + Cqs) (4-^32 — Fg\ — Fg2 — QG g3 — 2Fg3) 

~{Fg2 + G'g^ — 3Gg3 + G'g^){AFa2 — Fa2 + 36*2 — 66*03 ~ 2Fog3)} , 
a* = 2{Fg2 + G'g,-3Gg3 + G'g3) 

I \SFg1 + ~ 3G'c,3 + G'^3)(4Fo2 — Fa2 + 3(72 — 6Go3 — 2Fog3) 

— {Fca + 6*^4 — 36*03 + G'^.^[^Fg2 — Fgi — Fg2 — QGg^ — 2Fg3)^ , 
h = -\{{Fa2-Fa2W + {Fg2-Fg2-Fgi)g*^}. (3.6) 

To have a better feeling about the structure of the flow we plot the flow diagrams for i? ~ 1_ 
and -R ~ 1+ in the {g"^, a) plane (figure El ID) , based on the nature of the scaling dimensions we 
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discuss in the next sections. As we can extract all the quantities associated with the flow for 
any radius, in principle we can explore the behavior of the flow for a wide range of temperature. 
Remarkably we observe cyclic flow for R> 1.03, which is suggestive of resonances of high spin 
particles with a stringy behavior . The cyclic flow structure becomes very complicated, where 
the nontrivial flxed points show black hole like behavior. Note that a cyclic flow should give 
periodic c-function. However as the multi-critical points of matrix models are not necessarily 
unitary, this does not contradict with c-theorem. 




Figure 3: Schematic flow diagram for i? ~ 1_ in [g"^, a) plane with the flxed points Aq(0,0), 



3.2 The critical exponents 

Let us now go back to the matrix partition function we started with. After completing the RG 
transformations, it obeys the relation 



A*(0,a*), A*((7*2,0), A*((7*2,a*) 



Zn+i[9, a, R — 6R] ~ [A(A^, g,a, R — 6R\'^' Z^lg' = g + 6g, a = a + 6a, R] . 



(3.7) 



This leads to the Callan-Symanzik equation 



for the string partition function {i.e. the world-sheet free energy) 




r[g,a,R] 



(3.8) 



^[g, a, R] 



1 



\iaZ[g, a, R] . 



(3.9) 
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Figure 4: Schematic flow diagram for ~ 1+ in {g'^, a) plane with the flxed points Ao(0, 0), 
Al{0,a*), Al{g*^0), Al{g*^a*). 



The singular part of the world-sheet free energy is given by the solution of the homogeneous 
Callan-Symanzik equation. The inhomogenous part deflned by the change in the prefactor A, 
contributes to subtleties in the free energy. 

To discuss the critical exponents for the scahng variables, the renormalized bulk cosmological 
constant A = 1 — g/g*, and the renormalized fugacity for the vortices a — {1 — a /a*), we 
introduce the scaling dimension matrix 



k,l 



dAi 



(3.10) 



The eigenvalues Qj of this matrix represent the scaling dimensions of the relevant operators. 
The scaling dimensions at different fixed points are evaluated in the Appendix C. In terms of 
them, the homogeneous part of the Callan-Symanzik equation satisfied by the singular part of 
the free energy, around a fixed point, can be written as 



J^s[A,a,R] = 0. 



(3.11) 



The scahng of the free energy with respect to the renormalized cosmological constant, as one 
approaches the fixed point, goes as 

dR 



a 



InA-fii 



h{R)R 



(3.12) 
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where Fi is an arbitrary scaling function whose exphcit form depends on the initial conditions. 
Comparing the above expression of JF, with the matrix model result JF, ~ /\i'^~'ro) f^N'^/'y^ A], 
or using the standard definition of the susceptibility T ~ ^^\a=o ~ A^'''" , the string 
susceptibility exponent 70 is given by 



Note that in our analysis 2/71 ~ fii, i.e. 71 ~ 2/^7i is consistent with the matrix model relation 
7o + 7i = 2. This relation is independent of the explicit values of 70 and 71 and is easily 
obtainable from the consideration of the torus. The string susceptibility exponent at genus G 
is defined by 



3.3 Analysis of the flow 

Now we will analyze the nature of the fixed points considering the behavior of the scaling 
dimensions with respect to the change of the compactification radius R. Since the flow equations 
and the constraints together give rise to complicated system of nonlinear simultaneous equations, 
it would be easier to analyze them by studying the behavior of all quantities as functions 
of temperature. We will also plot the relevant functions for the convenience of the analysis. 
Referring to the Appendix C we observe that, the nontrivial fixed point produces a pair of 
c = 1 fixed points foT R — > 00, characterized by the asymptotic scaling dimensions {Qi = 1, ^2 = 
— 1.2} (more precisely, the asymptotic values are true for R > 1.08 and R> 2 respectively, see 
figure EI) and hence the string susceptibility exponent 70 = (2 — 2/Qi) = 0. 

These fixed points have one unstable direction. The critical coupling for such fixed points is 
vanishingly small {g*"^ 0,a* = 0) and hence the pair of c = 1 fixed points will be infinitesimally 
close to the gaussian fixed point, as we have already seen in the analysis of the ungauged model 
[T^ . As the radius is increased Qi grows to infinity as i? — » 0.70, and Q2 Aips sign from negative 
to positive at R = 1.03 and also grows to infinity as i? ^ 0.76. This indicates that the operator 
coupled to the fugacity of vortices becomes relevant at i? = 1.03 and triggers the expected 
liberation of the world-sheet vortices by BKT transition at the self dual radius [Til HSj. 

After this transition, the fixed point becomes purely repulsive up to -R ~ 0.67. In the range 
0.67 < R < 0.70, the critical coupling is infinitely large {g*'^ 00, a* = 0) as well as the scaling 
exponents Qi 00,^2 ^ 00. As the quantity h is also a large positive quantity here (figure 
Ej), the fixed point Al in this region is characterized by a negative specific heat, reminiscent of 
Euclidean black hole in fiat space-time. We will elaborate on this in the next section. These 
black hole like fixed points have a positive string susceptibility exponent 70 = (2 — 2/Qi) — 2. 



70 ~ (2 - 2/Q^) . 



(3.13) 



7g = 7o + G 71 . 



(3.14) 
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Figure 5: Behavior of the scahng dimension fin for with respect to x = irR 

Between < i? < 0.67 the matrix couphng g*"^, a* and both of the scahng dimensions 
are negative. Hence one gets a purely attractive fixed point with imaginary matrix couphng. 
However, as -R ^ 0, the couphngs are vanishingly smaU, also /i — 0, and one reaches a pure 
gravity fixed point (c = 0) with one unstable direction as Qi ~ 0.8 (such that 70 = 2—2.5 ~ —0.5) 
and Q2 —>■ —3.6. Note that even though Al exhibits so many important features for a wide range 
of temperature, the constraint we have here allows us to look at it only around R = 1 and around 
R = 0. We need to improve over this constraint. 

Similarly one can analyze the behavior of the fixed point A2 with respect to the parameter R. 
It exhibits black hole like behavior {h is large positive quantity) as < i? > 0.250. In this limit 
the fixed point A2{g* = 0,a* 00) is purely repulsive, Qi —*> 00,^2 00, and hence 70 — 2 
(see figure IE)). However, here the constraint on the radius is i? > 0.19. For the fixed points A3, 
there is also similar black hole like behavior with negative specific heat at 0.51 < R> 0.63 and 
at 0.12 < R> 0.19 and 70 = 2 (see figure El). There is no constraint on the radius. 

4 Comments On The Black Hole Fixed Point 

To study the thermodynamic behavior of the fixed points we would like to solve for the free 
energy from the Callan-Symanzik equation as one approaches the BKT phase transition (or 
the Hagedorn transition) and analyze its thermal properties. From the analysis of the previous 
section, one can see that the fixed point A* has rather exotic behavior around Rh ~ 0.73, little 



24 



Figure 6: Behavior of the scahng dimension 0,22 for with respect to x — nR 



befow (above) the Kosterhtz-Thouless radius (temperature) Rkt — l/(27rTxr) = 1- Around 
this region the scahng dimensions Qi, O2, and h, are large constants which simphfies the situa- 
tion. However, as we proceed we will see that the black hole like behavior would emerge from 
the region where /i is a large positive constant. Thus around Rh the free energy can be written 
as 



JF, ~ f[{R - RhIRhY"", a, a, N] , 



(4.1) 



where the inverse temperatures are defined by, 



(4.2) 



H 



The thermodynamic quantities are given by 



d{(5F) 
dp 



■PH)^hiZ{P-PH), 
-dlnZ 
dp ' 



E 



V 




(4.3) 



Using J^s, we have 
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Figure 7: Behavior of h for A.\ with respect to x = vri? 



E 



1 ^-1 



(4.4) 



Since near the phase transition the fluctuation of energy is large, the canonical ensemble would 
diverge. In such situation, it is better to pass to the microcanonical ensemble with fixed energy 
and the temperature defined by 



/3 



dS{E) 
dE 



Using ()4.4|1 for large positive h, one can solve for (3 in terms of E as 



E 



(4.5) 



(4.6) 



Combining this with the definition of temperature in the microcanonical ensemble ()4.5p one can 
calculate the near Hagedorn one loop finite energy correction to the usual definition of Hagedorn 
density of states, p{E) = exp[5'(£')] ~ explPnE], and the usual inverse temperature, which is 
otherwise a constant dS/dE = (3^. The finite energy corrections are of the form 
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Figure 8: Behavior of h for with respect to x = ttR 



dlnp ^ si 1 



O 



(4.7) 



Here the number si would come from the one loop correction. If si is negative, the specific heat 
is negative, i. e. increasing the energy of the system gives rise to the decrease of temperature, 
indicating Euclidean black hole like behavior in flat spacetime. 

Integrating our equation (j4.6|) to get the entropy and density of states we identify the one- 
loop correction as 

S{E) ~ PhE-^IhE, 
p{E) ~ E-Ti exp[/3jyE] , si = -l//i< . (4.8) 

This behavior is true for large positive h and the corresponding range of radius only. As we have 
analyzed in the previous section, we encounter such fixed points h.\ of very large and positive 
h in the region 0.67 < R < 0.70. These are pair of purely repulsive fixed points (over this 
region the scaling exponents Qi, Q2 are also large positive constants like h) of large (diverging) 
coupling g*"^ 00, a* = and positive string susceptibility exponent 70 ~ 2. We therefore 
identify such fixed points as the Euclidean Black hole like fixed points in the continuum limit. 
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Figure 9: Behavior of h for A3 with respect to x = ttR 



with a negative specific heat (figure ITn|l 

Using relations hke = — /?^f|j the entropy as a function of temperature is 

S{f3 - /?^) ~ -1 ( log(/5 - M - . (4.10) 

The discontinuity in entropy (a measure of the latent heat of the transition) suggest that the 
Hagedorn transition is a first order phase transition at little higher temperature than the KT 
temperature, driving the system to an unstable black hole phase. 



5 Discussions 

In this paper, starting from a matrix quantum mechanics on a circle with a periodic boundary 
condition (with gauge fields A{t) = Q''Q(t) as degrees of freedom apart from the matrix 
degrees of freedom 0(t) and also with the coupling a as the fugacity of vortices), we have 
analyzed the phase structure of the theory by a world-sheet renormalization group flow and 
observed following remarkable facts. The nontrivial fixed points of the flow does capture the 
known physics of the 2D string theory, namely a fixed point with the critical index of c = 1 for 
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Figure 10: Behavior of the specific heat (for Al) with respect to x = nR 



large radius (-R > 1), and moreover reveals a new phase. In our previous paper we have 
analyzed this c = 1 fixed point in detail showing that it exhibits expected logarithmic scaling 
violation of the singlet free energy and the T-duality. Note that for this class of fixed point (A^), 
a pair with a* = 0, the vanishing a forces A to be zero and the matrix quantum mechanics 
(jl.lll reduces to the usual matrix quantum mechanics with periodic boundary condition studied 
in The flow equations also exhibit indication of phase transition aX R = 1.03 as expected 

in BKT phase transition undergone by the 2D string theory due to liberation of world-sheet 
vortices. However, a new phase emerges beyond the phase transition at the self-dual radius. The 
fixed points Al exhibit an unstable black hole like phase with negative specific heat, between 
0.67 < R < 0.7. As -R ^ they end up into c = fixed points which is consistent with the 
expectation. For another class of fixed points A3, a pair with both nonzero a* and g*, there 
is unstable black hole like behavior with negative specific heat above the self-dual temperature 
(at 0.51 < R < 0.63 and at 0.12 < R < 0.19)^. From a thermodynamic study of the free energy 
obtained from the Callan-Symanzik equations we show that all these unstable phases do have 
negative specific heat. The thermodynamic quantities indicate that the system does undergo a 
finite temperature phase transition around the Hagedorn temperature {Rh = 0.7, around which 
the new phase is formed) and exhibits one loop finite energy correction to the Hagedorn density 
of states. 

Thus the thermodynamics of the 2D string theory above the BKT phase transition at the 
self-dual radius is governed by the unstable black hole like phase with negative specific heat 

^Working with Hl.l|l and then having a* = seems to improve the critical exponent in studying c = 1 fixed 
point with a nonzero h. However, we still need to improve on the constraint for A* which allows us to look at 
this particular fixed point around R = 1 and around R = only. 

■^Though for Aj, a pair of fixed points with nonzero a*, the behavior is seen for radii < i? < 0.25, there is 
a constraint to get the fixed point which renders R > 0.19. 
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for a range of temperatures beyond the BKT phase transition point. The remarkable thing is 
that we observe phases of negative specific heat for a range of temperatures rather than at a 
particular temperature like the case in ^Tj indicating that presumably there are many black 
holes at temperatures other than that of 2D cigar black hole described by dilaton gravity. It 
is then meaningful to ask where do the thermal black holes that we observe end up if they 
evaporate? Our analysis shows one class of these unstable objects evaporate to c = 0, the 
other classes end up in other black holes at higher temperatures. In 1 + 1 dimension, the string 
theory is integrable due to infinite number of conserved charges. Thus from continuum point 
of view it is possible that there are other black hole solutions not only characterized by mass 
and temperature but also by other values of conserved charges, in which case they might be at 
different temperatures. Here we can mention that the understanding that the integrability of 
the free fermion structure perhaps prevents the formation of c = 1 black hole is consistent from 
the nature of the flow and the fixed points. The c = 1 fixed points at large i?, dominated by 
the singlet sector, does not exhibit any black hole like behavior. 

To actually see that these objects are black holes and to deal with the questions like the 
formation and the unitarity of the scattering off the black hole one needs to study the dynamics 
of the nonsinglet states. As the observations indicate the nonsinglet states account for the 
entropy of the black hole, it would be nice to realize them as the excitations of the black 
hole. In black hole physics it was proposed that presumably the black hole entropy is due to 
open string with ends lying on the event horizon. It would be really interesting to study the 
nonsinglet boundaries of the matrix quantum mechanics in this context. The general framework 
of the world-sheet renormalization group approach is useful to study these objects, which is 
otherwise difficult. The Hamiltonian formalism, appropriate to address the questions of black 
hole dynamics in the Minkowski space, is complicated due to Calogero type of interaction. 
However, in the renormalization group approach one can utilize the Callan-Symanzic equations 
to calculate the wave functions for the nonsinglet boundaries or the macroscopic loops and can 
construct a ^-matrix. The wave equation for the macroscopic loops should then contain the 
information regarding the metric of the black hole to "see" the black hole at all. In this regard 
one nice observation from the renormalization group analysis is that the running of the radius 
with scale implies a deformation in the target space geometry if one considers the scale to go 
like dilaton. Presumably this could help to illuminate further the issue of getting the metric. 
Also it is interesting to understand the localized wave function (microscopic loop) describing 
the the tip of the cigar black hole |32! • The unique boundary associated to the inner core of the 
black hole is thought to be essential to understand the evaporation and the Hawking radiation 
j48j . We have work in progress on understanding the boundaries in this context. 

Regarding our RG method, as we have already discussed in jT^], it would be interesting to 
generalize the scheme for arbitrary couplings and to keep arbitrary powers of cj) in the evaluation 
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of the determinant obtained from the integration over the vector degrees of freedom {v, v*). Here 
we will just mention that (as discussed in our previous papaer) because of these vectors our 
partition function (in the a —* limit) essentially looks like the model discussed in 1^ 
and thus is useful to understand the presence of the boundaries. As a simple step towards 
generalization we would like to study the flow with an additional coupling due to 0^ term. 
These might reveal finer observations and also would be helpful to test the convergence of 
the scheme as well. The cyclic flow structure also deserves a detail study, specially near the 
regime of the black hole like behavior of the nontrivial fixed points where the flow structure 
is complicated. We would also like to understand how the relevant operators driving our flow 
look like in the matrix quantum mechanics side and how they translate to the operators in the 
world- sheet. 
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A The expression of a, ^at, An, R, N] 
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Evaluating the diagrams, T^lg, a, (f)N, A^, R, N] can be expressed as 
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B The Feynman Diagrams 



Here we evaluate and discuss the terms in different orders of the series a, ^at, A^v, -R, A^] 
()A.1|) using the summation rules discussed in section (3.3) and the relation ()2.10|) for the inverse 
Fourier Transform. 



B.l The terms of order 0($): 



g Tr 
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Now, changing the variables to 'center of mass' and 'relative' coordinates defined respectively 
by 

ti+t2 tl - t2 , , 
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we have, 
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and, 
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Hence, in 0{a), the coefficients of the 0($$) terms are given by comparing with the expression 
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where, after collecting all the results. 
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B.4 The terms of order O ($$$): 
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Using redefinition of tlie variables into tlie "center of mass" and tlie "relative coordinates", 
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where, 
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Now, 
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The last but one term can be evaluated as, 
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Now we will use the following redefinition of the variables into the "center of mass" and "relative 
coordinates" , 

T = ^{h + ^2 + h) , Ti = {h - t2) , T2 = {h - h) , 
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Now contribution of the different terms on the above sum can be evaluated as, 
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— exp ( — m(r2 — Ti)/R) ^exp ( — Uti/R^ exp (^imT2/R) — 2 exp ( — imri /i?) cxp(z/r2/i?))} 
+(ti + T2)| exp ( — in(T2 + Ti)/R) ^exp (ilri/R) exp (imT2/R) — 2 exp (irriTi/R) exp (ilT2/R)^ 

— exp (in(r2 + ti)/R) I exp ( — Uti/R) exp ( — imT2/R) — 2exp ( — imTi/ R) exp ( — ilT2/R 



a 



SnR 



exp{ilTi/ R) 
j^JpjR^ 



exp(— i/ri/i?) 



= — acothTri?. 
3 



(B.34) 



exp(i^(r2 -Ti)) 



/7ri? /'TtB, 
dri / (ir2(ri-T2) 
-TTil J-TVR 



E 



/Jtix / ■mT2. 



c^r2 J] 



/R^ 



exp ( — imT2/R) — exp ( — m(T2 — Ti) exp (i'mT2/R) exp ( — Uti/R) | 
+ (ti + r2)| exp {iWi/R) exp {imT2/R) exp ( - m(T2 + Ti)/i?) 

— exp ( — imT2/R) exp ( — Uti/R) exp (in(r2 + ri 
1 " " 



n - T2)| exp (m(r2 - Ti)/R) exp (Uti/R) 



E 



(iri 



exp ( - ilTi/R)^ + (ti + T2) (S{-Ti - T2)5'{t2) exp (Uti/R) 
-S{ti + T2)S'{-T2) exp ( - z/n/i?)) 



(n - r2) <5(r2 - n)6'{-r2) exp (^/n/i?) - 6{-T2 + ri)(5'(r2) 



(B.35) 
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{2.d) 



%a 



247r3i?3 



77/ 

exp(i-(T2-Ti)) 



(iri / o?T2(ti-T2) 



E 



exp {i—) exp ( - i^) 



R 



-ta 



247r3i?3 



dT2 



mn 



(ti — T2)| exp [ilTi/Rj exp ( — imT2/R) 



exp (m(r2 — Ti)/i?) — exp {iraT2/R) exp ( — i/ri/i?) exp ( — m(r2 — Ti)/i?) | 
+ + T2)| exp (ihi/R) exp (imT2/R) exp ( - m(T2 + Ti)/R) 

— exp ( — imT2/R) exp ( — Uti/R) exp (m(r2 + ri) /-R) | 

7"i - T2) ^5'(t2 - ri)(cosh(27ri? + r2 - ri) - cosh(r2 + Ti)) 



24sinh^ 7ri? Jo 
-5'{ti - T2)(cosh(27ri? + Ti - T2) - COSh(Ti + T2))j 

+ + r2) ^5'(-r2 - ri)(cosh(27ri? - T2 - n) - cosh(T2 - n)) - 

S'in + T2)(cosh(27ri? + Ti + r2) - C0Sh(T2 - Ti)) 

■ f cosh ttR H sinh 7ri?^ 

;V 6 12 / 



sinh 7ri? 



(B.36) 



(3) 



m.n.k 



(4n - 2m) /i? (2n - m)/i? | 



Hi _|_ 1 



rfti rft2 c^ts Tr(A(ti)0(t2)0(t3)) Yl 



-^m—n—k4'k4'n—m 

{4n-2m)/R {2n-m)/R 



m,n,K 

ex.p{ik{ti — t2)/R) exp{im{t3 — ti)/R) exp(in(ti — t^)/R) 

r2irR rirR 



ml _|_ 1 



ni _|_ 1 



'0 J--KR 

exp{ik{T2 — Ti) / R) exp{imTi / R) exp(— inri/i?) 
= 0. 



(B.37) 



As the overall behavior of the function is proportional to dri dT2 (ri — T2)S{ti — T2) , 

therefore the contribution vanishes. Similarly the contribution of the other term with similar 
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sum 



^« \ V~ ~ . r ^m-n-k<t>n-m<t>k (B.38) 



— + 1 — + 1 f 



vanishes also. Thus, the 0{A[(f), 0]) term is given by 

(G,3 / + G'c.s «) j dt A{t)[ct>{t), ^{t)] , 



where, 



Gg3{R) = ^(7ri?cosh27ri? - ^ sinh27ri?) 

J. ^ ^ 



G^gfi?) = --cothnR \ (— cosh7ri?+ — sinhTTi?) . (B.39) 

3 sinh Trit! 6 12 



B.7 The terms of order 0(A2$2) 



(1). 2a Tr 



[n 



2nm + Aral - 2nl)/R^ 



m,n,l,k (^2 + l) (ij2 + l) 

f27ri? 



^'Pl-n'Pn-mA.m-l—kA.k 



2a f 

——4 / dti dt2 dh dU TT{A{ti)<P{t2)<i>ih)A{U)) V 



(n^ - 2nm + 4ml - 2nl) /R^ 



Aj / T\i TJX 

exp(i— (ts - ^4)) exp(i— (ta -ti)) exp(i— (ti - ^2)) exp(i— (tg - ^3)) 



(B.40) 



Again using the usual redefinition of the variables into the "center of mass" and "relative 
coordinates" , 



r= -(^1 + ^2 + ^3 + ^4), Ti = (i4-ti), T2^{u-t2) n^iu-ta), 

dti dt2 dt^ dti — j( 2,' ^' ^' dT dri dT2 dr^ — dT dri dr2 dr^ . 

\T,Ti,T2,r^/ 

Considering ri, T2 and T4, to be small and keeping terms of the form 0(y4(t)^$(t)^) only, above 
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series could be evaluated as, 



2a Tr 

2a 



E 



(n^ - 2nm + Ami - 2nl)/R^ 



{27r Ry 



/ttR 
dri dT2 ^ 
"^^ m,n,l 



iin?-2nm + Aml-2nl)/R^ 



cxp(!-^(r2-Ti)) exp(-i^(T2)) / dT Tv{,ti'-iT)A\T) 



2nR 



4a 



2-kR 



dT T^((t>\T)A^{T)) 



(n^ - 2nm + Ami - 2nl) jR? 



(27ri?)3 Jo 

'^^ I n m 

dn dT2 exp(i— (fi + f2)) exp(-i— 2f2) exp(i — (f2 - n)) 



:b.4i) 



The terms containing powers of n in the numerator of the sum inserts higher and higher deriva- 
tives of delta function over time variables and are eventually computed to be zero. The only 
non- vanishing contribution comes from the term without n in the numerator. 



4 a 
(27rRy 



2nR 



r f I n m 

J dri dr2 exp(i— (fi + fs)) exp(-ii— 2f2) exp(i— (fa - n)) 

I 77/ Tft I 

+ exp(i— (fi - n)) exp(i— (2f2)) exp(-i— (f2 + fi)) + exp(-i— (fi - T2)) 



Tl Tfl I Ti Tft \ 

exp(i— (-2f2)) exp(i— (f2 + fi)) + exp(-i— (fi + f2)) exp(i— (2x2) exp(-i — (f2 - fi))j 
j dT Tt(^(P\T)A^{T)^ J dfi(^sinh(7ri?-fi)sinh(7ri? + fi)) 

1 r^""^ / \ 

;(7ri?cosh27ri?- -sinh27ri?) / dT TrU'^{T)A\T)) (B.42) 



16 a 
sinh^ nR 

8 a 



sinh ttR 
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(2) Similarly, 



-2 a Tr 



E {2171^ — nm + 2ml — nl) / . 
f „,2 , -, \ / p — {(Pn-m^m-l^l-n-k<Pk) 



m,n,l,k 
■27rR 



2 a f 

- / ^^3 dU Tr{(f){ti)A{t2)A{ts)(f){U)) 

exp(i— (^3 - ^4)) exp(i— (t2 - ^3)) exp(i— (^3 - ti)) exp(i — (ti - ta)) 



(2m^ — nm + 2m/ — nl)/ B? 



2 a 



(27r7?) 



-ttR 



m,n,l,k 



exp(-i— Ts) exp(i— n) exp(i— (ts - ti)) 



~ a(-cotli7ri2-27ri?(l + cotliVi?)) / dT Tr{(j){T)A{T)A{T)(j){T)) (B.43) 

In evaluating above expression again we see that the term containing (any power of) n in the 
numerator of the sum is not contributing. Again, 



2 q; Tr 



{2f-nm + 2ml-nl)/R^ 

/ ^ 77;;! -,\ r 1.2 {Am-lVl-nVkAn-m-k) 



. m,n,l,k 

2 a '•^-^ 



(27ri?) ./o 



, / < / N , / N , / N . / N \ (2/^ — nm + 2ml — nl)/R 

dtidt2dhdUTl{A{ti)(P{t2)(p{h)A{U)) 22 7^2 TTTp 



m,n,l,k 



k I 7X TYt 

exp(i— (^4 - ts)) exp(i— (ti -ta)) exp(i— (ts - ^4)) exp(i— (^4 - ti)) 



— — / dT Tr{A{T)<f>{T)<l>{T)A{T)) / dndr^ J] 



exp(i— Ti) exp(-i— T2) exp(i— (t2 - ti)) . 



(2/2 - nm + 2m/ - nl) / R^ 



(B.44) 



Exchanging / -^-^ m and ri ^ —T2, this gives equal and opposite contribution to the previous 
expression and hence all the similar pair of sums 



(2m2 -nm + 2ml - nl) /R^ 

m,n,l,k + + 

(2/2 - /,/;/ + 2//// - n1)/R^ 



E 

n,n,l, 



{2Ai_n-k<Pk4>n-mAm-l — 4>kAl- 



■n—k'rn—m 



Am-l) 



m,n,l,k ('r2 + 1)(r2 + 1) 

give zero contribution. 



{'^4'l-n4'kAn-m-kAm-l ~ 4>l-nAn-m-k'PkAm-l) i 

(B.45) 
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(3) a Tr 



E 



4 1 

+ 



a 



m,n,k,k' ^ R 
2-kR 



+ 1 



+ 1 



m—k' 



(27ri?)4 Jo 



dii dt2 dts dU T^{A{ti)(t){t2)(t>{h)A{U)) J] 



m,n,k,k' R"^ 



i?2 



exp(i— (ti - ts)) exp(i— (^4 - t^)) ex.p(i—{ti - ^4)) exp(i— (t4 - ti)) 



~ — cothTTi? / dr Tr(A(r)(/)(T)(/)(r)A(T)) 
2 Jo 

(4) Similraly, 



(B.46) 



—a Tr 



— + 1 

m,n,k,k' R^ 

^2-kR 



{4'kAm-n-k4'k' An-m-k' + 4'kAn-m-k'4'k' A^-n-k) 



= / dU Tr{A{h)<P{t2)A{t3)(l){U)) (^TT 

|exp(i— (i2 - ii)) exp(i— (^4 -ig)) exp(i— (^2 - ^4)) exp(i— (^4 - ^2)) 



+ exp(i— (t4 - ti)) exp(i— (t2 - ta)) exp(i— (t4 - t2)) exp(i— (ta - t4))| 

/■27ril 

~ -2Q;coth7ri? / dT Tr(A(r)0(r)A(r)0(r)) 
Jo 



(B.47) 



(5) 



2(27ri?)4 



-n0n— m^m— Z— fe^fe j 



/■27ril 

/ dti dt2 dh dU Ti{A{h)(t){t2)A{h)(t){U)) 



m.n.l.k 



exp(i— (^3 - i4)) exp(i— (ta - ii)) exp(i— (ii - ^2)) exp(i— (^2 - ^3)) 



16(7ri?)- ^0 



27rR 



dnk(A(T)MT)A(T)MT)) 



I dridr2 exp(-i— r2) exp(i— ri) exp(i — (t2 - Ti)) , 

7-7rR K K K 



(B.48) 
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where, 



E 



tB 



I Ti in 

NO / / 9 X / 9 X / ,9 X I dTidT2 exp(— i— T2) expfi— Ti) expfi — (t2 — ti)) 



^7^ 



87ri?sinh vri? 



X/ TT;:^ 7Y / ^'^i^'^2 ( sinh(7ri? + T2) smh(7ri? - n) exp(i — (ts - Ti)) + 

m [w + ^) Jo ^ 



smh(7ri? + T2) sinh(7ri? + ti) exp(i— (t2 + Ti))) 



167ri? sinhTri? 

Tfl 

exp(i— (t-2 - ri)) + 



T^^^a / dTxdr^ ( ( cosh(27ri? + r2 - n) - cosh(ri + r2)) 

^ + ij ^0 ^ 



( cosh(27ri? + Ti + r2) - cosh(Ti - T2)) exp(i— (r2 + Ti)) j 

V -r^ / dT df [{ cosh(27ri? + 2f) - cosh(2r)) exp(i— f) + 



8nR sinh 7ri? 



.2m, 



(cosh(27ri? + 2T) - cosh(2f)) exp(«— T) 

{ttR cosh(27ri? + 2f) - ^ sinh(27ri?)) J] . , ' x exp(i^f) + 



SttR sinh vri? 



2m 



■kR/2 



ttR 



J dT {nRcosh{2nR + 2T) - sinh(7ri?)) ^.^j^^^ 



cosh(7ri? - 2T) 



9 



8 sinh^ 71 R 



dT 







,7rR 



(— cosh(27ri? + - sinh 27ri?) cosh(7ri? - T) + (— cosh(27ri? - T) 

^4 ^ 



-- sinh 27ri?) cosh(7ri? + T) + {■KRcosh{2nR + 2T) - sinh(7ri?)) cosh(7ri? - 2T) 



9' 



1 1 

TT R cosh 37rii H sinh 57rii H sinh 37rii 

8 4 



8sinh^ ttR 

3 X 
— sinh ttR + - sinh^ 27ri? — sinh^ ttR ] 
8 4 / 



(B.49) 



Collecting all the terms, the 0{A(j)A(f)) and the 0{(f)'^A'^) terms are given by 

iGg49' + Go^^a J dt TvA{t)ct>{t)A{t)ct>{t) + {G'^^g^ + G:,4a) ^ dt TvA\t)ct>''{t) , 
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where, 



GndiR) — \ (tt^R^ cosh SttR + sinh 5t:R + sinh SttR — - sinh ttR 

^ Ssinh^TTi^V 8 4 8 

+- sinh^ 27: R - sinh^ ttR] 
4 / 

G'g,iR) = 

Gad{R) = -2coth7ri? 

G'^^{R) = - coth ttR . (B.50) 

C The Scaling Dimensions 

For A*, 

_ 4Fg3 + 3Fg2 + 3Fgi + ^Ggi 
"11 — ^ ) 

Fg2 — Fgl — Fg2 — 5Gg4 — 2Fg3 



3Fn2 + 2Fn 



12 



a2 ^-rag3 



Fg2 — Fgl — Fg2 — 3Gg4 " 2 Fg^ 

^21^0, 



22 



^ _2 / -^g2 + Fg2 + Fgl - 2Gg3 + |Gg4 + Fg3 \ 

^ Fg2 — -Fgl — Fg2 — 3Gg4 " 2 Fg^ ^ 



For A^, 



_ 6Fq,(,3 + 6Fa2 + 860,4 — 2G2 

Sill — ~ ) 

Foa — Fa2 — 3Ga4 + 2G2 — 2Fag3 

^12^0, 
^21^0, 

^ _ 6Fq2 + 2Fa2 — SGgS + 36*04 + 2G2 + 2FQ,g3 
-Fq2 — Fa2 — 36^4 + 2G2 — "^FagS 
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